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The relation between the interaction parameters for fermions on the spatial lattice and 
the two-body T matrix is discussed. The presented method allows determination of 
the interaction parameters through the relatively simple computational scheme which 
include the effect of finite lattice spacing. In particular the relation between the inter- 
^ ' action parameters and the effective range expansion parameters is derived in the limit 

' of large lattices. 

in ■ 

1. Introduction 

One of the most promising approaches to describe the nonrelativistic many-fermion 
system is to define the problem on the spatial lattice and subsequently apply the 
. Monte Carlo technique to probe the configuration space. This approach is particu- 

larly efficient when dealing with strongly interacting and dilute systems. In the limit 
, , of low densities the interaction between particles has predominantly two-body char- 

rS , acter. Moreover it can be described by few parameters only, characterizing the low 

' energy physics of two-body collisions. For example in the case of trapped fermionic 

atoms, which have been recently the subject of intensive theoretical studies (see^ 
and references therein), the scattering length is the only parameter which deter- 
mines the interatomic interaction. Similarly, the physics of dilute neutron matter at 
densities corresponding to kp < 0.6 fm^^ is captured by two parameters: scattering 
length and effective range in the ^5*0 channel. 

The basic problem of the lattice formulation is the determination of the Hamilto- 
nian parameters in order to describe properly the scattering phase shifts of two-body 
collisions. The prescription given by Liischei^l relates the discrete energy spectrum 
of two-particle states in a box to the scattering phase shifts. However in order to 
use the Liischer's formula, the exact two-particle energy spectrum has to be deter- 
mined first, which in general is a nontrivial task. Another prescription, given by 
Seki and Kolck for the case of large lattices, is more suitable for the low-energy 
nuclear problems. In this approach the interaction parameters are determined by 
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consistently applying the effective field theory power counting rules. 

In this paper we present an alternative approach to determine the two-body 
interaction parameters. The method is especially well suited to deal with the dis- 
crete form of the interaction, commonly used in the lattice calculations .l^fSiSi p^\\ 
finite lattice spacing effects are properly included in the limit of large lattice vol- 
ume. The final prescription bypass the tedious Liischer's algebra leading to a simple 
expression convenient for practical applications. 

2. Interaction on the lattice 

To capture the low energy physics of a Fermi system let us consider the interaction 
which depends only on the relative coordinates of two particles: 

V = \Y. /rf'rdVV'l(r)7^1,(r')T/(r-r')V'A'(r')^A(r), (1) 

where the field operators obey the fermionic anticommutation relations 
{'0|(r), -0^' (r')} = '5aa"5(i" — r') and A denotes the spin degree of freedom. In order 
to place the system on a 3D spatial lattice with lattice spacing h and size L — Ngb 
the discretization of the above expression has to be performed. This leads to the 
replacement: ip\{r) iip\{ri), where ri = ib and i = {ix,iy,iz) denotes lattice sites, 
and ix,iy,iz = 1, . . . ,Ns. Consequently one gets instead of ([T]): 

^^l^' E E^l(^i)^l'("'j)^('-i-^j)'^>'('-j)^A(ri). (2) 

A,A' = Ti ij 

Note that for the problem defined on the lattice, we need only a discrete set of 
values of the interaction potential defined for the lattice sites Rfe(i,j) = ri — rj. 
Hence it is sufficient to introduce the interaction on the lattice which is defined 
only on the lattice sites. It can be written in the form: 

D 

V{r-r') = Y,9kS{r-r' -Ilk). (3) 

fe=i 

The coupling constants gk contain the full information about the interaction be- 
tween the particles. The simplest case corresponds to the on-site interaction when g^ 
is nonzero only when = {D = 1). A more complicated case when the neighbor- 
ing sites are included leads to the so called extended Hubbard model {D = 7), widely 
used in Monte Carlo simulations (see^and references therein). It requires effectively 
two coupling constants (due to the rotational invariance of the Hamiltonian) : go for 
R = (0,0,0) (on-site interaction) and gi for R e {(±6, 0, 0), (0, ±6, 0), (0, 0, ±6)} 
(nearest- neighbor interaction). Then the interaction reads: 

V^gob'J2^^{rOHin) + ^b' E ^^^^('•^^^'(^j)' 

i A,A' = Ti <ij> 

where n\{ri) = ip\{ri)ijj\{ri), and < i,j > denotes the summation over all neigh- 
boring pairs. 



I 



November 7, 2009 20:19 WSPC/INSTRUCTION FILE 

Proceeding ' Kazimierz '2009 



Effective range expansion for the interaction defined on the lattice 3 

3. T matrix for the lattice interaction 

In order to link the coupling constants gt to the phase shifts one has to consider 
the T matrix, which in the case of low energy scattering can be parametrized as: 

- — 7p-p ~~- + ^ W - *P + 0(p4), (5) 
m a 2 

where a is the scattering length and rcff denotes the effective range. The T matrix 
fulfills the Lippmann-Schwinger equation: 

Tpp' = Vpp' + -jj^ V^pkGp'kTkp' , (6) 

k 

where Gpp' = m/{p'^ — p'^) is the free particle propagator with the reduced 
mass m/2. Summation over momenta is limited to the first Brillouin zone ie. 
—n/b ^ Px,y,z ^ T^/b, and each momentum component is discretized due to the box 
quantization conditions. As a first step, let us rewrite Eq. ^ in an iterative form, 
more suitable for numerical applications. Expanding the Lippmann-Schwinger equa- 
tion one obtains: 

J"pp' = Vpp' + jj ^ VpkGp'kVkp' + jj ^ ^pkGp'kVkk'Gp/k' Vk'p/ -I- . . . (7) 

k k,k' 



If we define the matrix M*^^) as: 

3- ^ E ^pi^Gp'kAC (8) 



L3 

k 



where = Vpp' then Eq. ([T]) takes the form: 



pp 



Tpp' = M^J, + M« ^pkGp'kA^ + • ■ • (9) 

k 

Clearly the T matrix can be written as an infinite sum: 

oc 

7^PP'=E<P" (10) 
where the M*^") matrices are related to each other through the recurrence relation: 

Kl''^-]iT.Vp^Gp,^<i'^ (11) 

k 

Kp'=^pp'- (12) 
The equation Eq. PH)) holds in general case. However in the case when the effective 
range expansion Eq. ([5]) is valid, only the diagonal matrix elements are required. For 
the lattice interaction given by Eq. ([3]) the contribution M*^°) takes the following 
form: 

D 

^4p = E 5«e-*P^'e'P^- = Tr L^°^ , (13) 

i=l 
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where the square matrix L'-*'-' of dimension D x D is defined as: 



L(0)=5,e-^P(^-^^^ (14) 



The contribution to M^^^ reads 

D 



^4p^ - E .9.5.e-P(^'-^^)Gp(R, -R,) 

D 

= E ajLfG^iR, R.) = TrL(i), (15) 



where: 



Li]' =^5fei(°)Gp(R,-R,), (16) 

GpW = i3E^"''^""^pk- (17) 

k 

Continuing this procedure for higher order contributions it can be shown that: 

oo 

Tpp^Y.^TL^"\ (18) 

71=0 

where: 



4"^'^=E.9^4'^Gp(R.-R,), (19) 



(n+l) - ^ r(")^ 
fc=l 

Lg) =5,e-^P(^-^^). (20) 



Hence the problem of computing diagonal elements of the T matrix was reduced to 
calculation of traces of a relatively small matrices L^'^\ of dimension D x D. It is 
important to note that the final result (|18p includes the effects related to the finite 
lattice spacing (incorporated by the relative coordinates of the lattice sites Rfc). 

Let us examine the large lattice limit, by letting L — > oo while keeping the lat- 
tice constant b fixed. Then the momentum is continuous within the first Brillouin 
zone, and the summation can be replaced by the integration. To simplify the anal- 
ysis, however, we place a spherically symmetric cut-off, including only momenta 
satisfying p < pcut = T^/b: 

tfp dp^ dpe cospe / dpp^. (21) 

Jo Jo Jo 

This prescription sets to zero all two-body matrix elements, if the relative momen- 
tum of two particles exceeds a given momentum cut-off. Note that the function 
Gp(r) will depend on the momentum cut-off and consequently the scattering pa- 
rameters like a i Toff will also be the functions of Pcut- Performing the integration 
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with respect to the variables dp^ i dpg the Eq. p7p transforms into: 

m Z'^''"' sinfcr 
/ dkk— 



G^ir)^^l dkk-^. (22) 



The remaining integral can be calculated analytically by expanding sin kr and using 
the relation: 

-^^ = P--tnS{x), (23) 
X + iO+ X 

where V stands for the principal value. The final result has the form: 

G^i^) = -^j:y-Uy,F(j^p) (24) 



where: 



1 n 21 2j-2l + l _ 
^ 2j - 2/ + 1 Pent + P 

The expression (|18p together with ([24]) and ([5]) provides the most convenient pre- 
scription for numerical applications. 

4. Example: on-site interaction 

As an example let us consider the on-site interaction: 

V{r - r') = god{r - r'). (26) 



It corresponds to the well-known Hubbard model ar id is p resently widely used to 
simulate the system of dilute, cold fermionic atoms! I I I In this particular case, 
it is possible to find an analytic formula relating the coupling constant go to the 
scattering length and the effective range. It is easy to realize that L*^"'' become now 
one-dimensional matrices ie.: 

i^"^ = 30 (ffoGp)", (27) 

where: 

2pcut - P In I 1 

47r^ \ p- Pent 

The sum of the geometric series reads: 



Gp = Gp(0) = -— ( 2pcut - P In le-^l + ) • (28) 



In order to find the inverse of the diagonal matrix elements of T for small values of 
the momentum, we expand the logarithmic function: 

In « (30) 



I 



I 



November 7, 2009 20:19 WSPC/INSTRUCTION FILE 

Proceeding ' Kazimierz '2009 



6 G. Wlazlowski, P. Magierski 

and get: 

T = \ p -tp+0{p). 31) 

m mgo tt Trpcut 

Comparing this equation with Eq. ([5|) one reproduces the weU-known resultJ^ 
1 47r 2pcut 4 

- = h , rcff= . (32) 

a mgo n Trpcut 

It is worth noting that due to the finite lattice spacing the zero range potential 
acquires the non-zero effective range red- 



5. Conclusions 

In this work we have presented the method to determine the two-body effective 
interaction parameters from the low energy scattering data when the problem is 
defined on the lattice. In the limit of large lattices we have obtained expressions 
which are particularly convenient for numerical applications. The prescription re- 
lates the interaction coupling constants to the effective range expansion parameters 
and can be applied to Monte Carlo simulations of many-body systems on large 
lattices. 
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